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HIGHLIGHTS 





è An aging based optimal scheduling framework is developed. 

è Operating conditions and maintenance intervals are optimized simultaneously. 

è Aging cost is considered in the objective function. 

è Equivalent operating hour approach is proposed to address component aging effects. 
è Results show 4% higher annual profit for a gas turbine power plant. 





ARTICLE INFO ABSTRACT 





In this paper a scheduling optimization framework is developed to enhance power plants operational decision 
making process. The proposed framework optimizes plant schedule including operating conditions and main- 
tenance intervals simultaneously and on an hourly basis. In a long term operation plant performance deteriorates 
due to components aging. This study employs equivalent operating hour (EOH) approach to describe compo- 
nents aging impact on the plant performance deterioration and consequently plant long term profit. Modeling of 
components aging increases system simulation accuracy in long term operation and the optimum decision 
variables would be more reliable and realistic. 

Validity and usefulness of the proposed methodology are demonstrated by optimizing the operating condi- 
tions and maintenance intervals of a gas turbine power plant, under different seasonal ambient conditions and 
energy prices. The case study results effectively meet all the positive expectation that are placed on the proposed 
aging based optimal scheduling framework. Results show that optimal operation schedule depends on the 
maintenance intervals. Therefore, operating conditions and maintenance intervals should be optimized si- 
multaneously. In addition, derived optimal schedule increases the system’s long term profit approximately 4% 
annually in comparison with a case that operates at optimal schedule without considering the aging effects in the 
optimization procedure. 
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1. Introduction expected total cost of plant due to operation, failures and preventive 
maintenances. 
In recent years, the use of optimization models to determine plant 


optimal schedule has earned popularity as in [1]. Scheduling is widely 


1.1. Motivation 


Besides using high-tech components in power plants, plant opera- 
tion optimization can significantly improve energy efficiency and eco- 
nomic performance, as efficiency of plant components generally de- 
pends on operating conditions. In addition, system preventive 
maintenance can reduce plant operation and failure costs, however it is 
also costly when done frequently. Therefore, optimizing operating 
conditions and preventive maintenance intervals can minimize the 
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used to maintain and establish operating conditions and maintenance 
intervals of a plant over time. It should be noted that power plant 
components are degraded through long term operation [2]. Therefore, 
components performance profile over time varies. To have a reliable 
and accurate scheduling optimization results it is necessary to consider 
components aging models in the optimization procedure. The effect of 
considering degradation mechanisms in the optimization procedure is 
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Nomenclature 

a temperature correction factor of p 
ar’ temperature correction factor of hr 
b steam injection correction factor of p 
E steam injection correction factor of hr 
6 load correction factor of p 

ô’ load correction factor of hr 

n efficiency 

A aging correction factor of p 

NV aging correction factor of hr 

A ambient condition 

G cost 

Ce levelized aging cost 

c fuel cost 

C heat price 

CG maintenance cost 

CG electricity price 

d day number 

dpe performance deterioration 

EOH equivalent operating hour 

EOH failure EOH 

F fuel consumption 

h useful output heat 

hr heat rate 

HR? heat rate standard value 

I model input 
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kok water/steam injection severity factor 

mf maintenance factor 

mf" maintenance factor of combustion inspection 
mf maintenance factor of hot gas path inspection 
mf maintenance factor of rotor inspection 

Ms steam percentage to compressor inlet air flow 
O model output 

op operating condition 

op! operating condition lower band 

op" operating condition upper band 

p output power 

P system performance 

p° base load output power 

R income 

S load severity factor 

S operating schedule 

t hours of system operation 

T temperature 

g long term profit 

Subscripts 

i hour 

j month 

k year 

m maintenance type 








investigated in [3] and the degradation based optimization (DBO) 
concept is introduced [3]. The combination of plant optimal scheduling 
and aging models is an extended approach that is presented in our study 
and a framework is developed. The proposed framework optimizes 
plant schedule over long term horizon considering components aging. 
The framework outputs determine the plant startup time, production 
level and maintenance intervals. This framework can be used in the 
sensitivity analysis of energy price and ambient conditions as well (see 
Fig. 1). 


1.2. Literature review 


Several current research address power plant scheduling without 
considering plant aging [4-6]. In these studies, mathematical ap- 
proaches for optimizing power plant schedule according to different 
ambient conditions (temperature, humidity, ...), system characteristic 
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and economic assumptions (energy and fuel prices, and fluctuating 
energy sale prices) are available. 

Other studies have developed aging models that estimate long-term 
power plant aging and performance deterioration. Aging models of 
power plant components is a well-developed field [7-9]. However, little 
emphasis has been placed on combining operation scheduling optimi- 
zation and plant aging models. 

Related pioneering works which introduce plants performance de- 
terioration into scheduling, are described in [10,11]. These references 
focus on considering the effects of operation conditions on system 
maintenance interval and failure rate. In [12], the intention is to 
compute power plant equivalent operating hours (EOH) that determine 
component aging and performance deterioration under standard oper- 
ating conditions. This concept is very appealing, however, specific dy- 
namic lifetime scheduling model is not developed. In our study, the 
framework for a hybrid approach, considering performance 
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Fig. 1. Data flow schematic of aging based optimal scheduling. 


T. Parhizkar et al. 


deterioration in operation scheduling is developed to have a large im- 
pact on a power plant’s economic performance. 

The proposed framework, is utilized to cope with the problem of 
hourly operational management of power plants that exchange elec- 
trical energy with the grid, whose prices change hourly and day by day. 
Hourly long term optimization makes the optimization procedure more 
complex. Several optimization algorithms of scheduling can be found in 
literature [13,14]. In [15] linear model is developed and hourly op- 
timum operating conditions are derived based on two stage optimiza- 
tion procedure. Time horizon at the first stage is hourly and then at the 
second stage is monthly. Many of the studies use linear approximations 
to describe inherently non-linear energy conversion processes [16,17], 
while others overcome this approximation through non-linear and 
mixed integer non-linear programming [18,6,19,20]. 

As can be concluded from the reviewed literature, most of the re- 
viewed literatures did not consider the component aging effects in the 
scheduling optimization procedure. In addition, most of the studies 
focus on the operation scheduling or maintenance intervals optimiza- 
tion, while none of them has considered the interaction of these two 
aspects simultaneously. 
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1.3. Paper contributions 


The major novelties of this study include the following: 


One of the main novelties of the proposed framework is the con- 
sideration of aging effects in the optimization procedure. The 
components aging can result in performance deterioration and 
system failure. The rate of aging process depends on many factors 
including operating conditions, system characteristics and period of 
operation. Effective management of components aging is a key 
element of the efficient and reliable operation of plants. The de- 
veloped framework manages aging process by optimizing system’s 
schedule in order to maximize the lifetime profit. 

An improved EOH approach is used in order to address components 
aging effect on plant performance deterioration. This new approach 
results in more accuracy in long term simulation and more reliable 
and realistic plant optimal schedule. 

The developed aging based optimal scheduling framework considers 
aging cost in the objective function as well as components aging in 
the optimization procedure. Aging cost is defined as the hourly 
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Fig. 2. Integrated framework of aging based optimal scheduling. 
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Fig. 3. An I/O model flow data. 


preventive and corrective maintenance costs. As a result, optimal 
hourly schedule is affected by not only the income of selling elec- 
tricity and operation cost but also the maintenance cost. Therefore, 
plant hourly profit is more realistic and the optimal schedule has a 
higher lifetime profit in comparison with other scheduling meth- 
odologies such as day-ahead method. 

e The objective of the proposed framework is to find system optimal 
schedule on an hourly basis in long term operation that is a new 
approach in schedule optimization. The model inputs such as am- 
bient conditions, fuel and electricity price varies on an hourly basis 
and consequently optimal schedule will change hourly in long term. 

e The proposed framework produces optimal maintenance and op- 
eration intervals of power plants simultaneously. As maintenance 
intervals and operation schedule are interrelated, the model opti- 
mizes both decision variables simultaneously and optimum opera- 
tion and maintenance intervals schedule will be derived at once. 

e In order to optimize plant hourly schedule and maintenance inter- 
vals, an innovative two stage optimization algorithm is proposed. 
The proposed heuristic algorithm is the hybrid of gradient search 
method and branch-and-reduce optimization navigator (BARON) 
method. The low computation time of this method allow concluding 
that this new optimization algorithm has a high performance in 
scheduling problems. 

e In addition, a 9F.05 model of general electric (GE) gas turbine 
power plant is modeled in this study. The hourly optimal schedule 
and maintenance intervals over 15 years of operation is derived for 
this power plant. Results show the effectiveness of the proposed 
framework and solution algorithm. 


1.4. Paper organization 


The paper is organized as follows: A brief overview of the developed 
framework is introduced in Section 2. The framework modeling and 
formulation are presented in Section 3. In Section 4 the optimization 
algorithm is fully described. Section 5 presents the characteristics and 
model of the system under study. In Section 6, the results of the case 
study are discussed, and finally conclusions and contributions are 
drawn in Section 7. 
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Fig. 5. BARON method algorithm. 


2. Aging based optimal scheduling framework 
2.1. General description 


The proposed approach aims at deriving maximum profit schedule 
that meet power plant operational and economic constraints. This ap- 
proach considers the plant performance deterioration due to aging ef- 
fects in the optimization procedure. Therefore, an aging based optimi- 
zation framework is proposed, that offers a comprehensive 
representation of the main factors of aging processes in power plants. 

Aging based scheduling model is categorized as a dynamic non- 
linear programming model. Such a model consists of three main ele- 
ments: (1) an objective function, expressing the power plant lifetime 
profit to be maximized; (2) constraints, equality and inequality equa- 
tions expressing the power plant limitations and logical relationships 
that must be satisfied; and (3) decision variables, the unknown schedule 
to be determined by the optimization. 


Fig. 4. Flow data of I/O and aging models integration. 
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Fig. 6. Optimization algorithm of aging based 
optimal scheduling. 


] 
] 
] 
] 
] 
i 
i 
] 
l 
-l 


——> Deteriorated performance 













No Satisfy the 


stopping 
criterion? 


Yes 





The objective function of the model is to maximize power plant 
lifetime profit. Profit is defined as income from selling power plant 
products such as electricity and heat, minus system total costs including 
fuel and aging costs. The aging cost is the hourly maintenance cost of 
the system as a function of plant schedule. 

Several constraints must be satisfied in the optimization procedure. 
The equality constraints are mostly the technical governing equations 
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of the power plant components and processes. The inequality con- 
straints are the economical and operational limits of the plant. 

The decision variable of the model is system schedule that includes 
system operating conditions and outage plan of the power plant. Other 
main model outputs are plant hourly profit, income and costs including 
fuel and maintenance cost, as well as hourly and accumulated plant 
performance deterioration. 
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Air 


| 






Combustion 
chamber 


Exhaust temperature 640 °C 
Exhaust energy 1794 MM kJ/hr 


Table 1 
Technical characteristics of the studied gas turbine. 


Manufacturer Model Power output Heat rate Nominalspeed Cost 
= z MW kJ/kWh rpm $/kW 
GE 9F.05 299 9295 3600 183.14 


Finally, the modified BARON global optimization algorithm is pro- 
posed to derive optimal plant schedule that maximize the objective 
function while respecting all the model constraints. The solution 
methodology was explained in more detail in the revised manuscript as 
well. 


2.2. Flaw diagram 


The flaw diagram of the aging based optimal scheduling framework 
is shown in Fig. 2. The framework determines optimum hourly opera- 
tion and maintenance schedules based on the hourly weather, hourly 
energy prices, and operational/economic constraints. 

The optimization model is the main part of the framework that 
consist of two main models. The components input/output (1/0) model 
focuses on system performance and efficiency. The component aging 
model considers the updates of system characteristic due to degradation 
mechanisms which is correspond to the operating conditions. 

To consider the effect of operating condition variations on perfor- 
mance deterioration, the maintenance factor and EOH concepts are 
introduced [12]. conditions and full load operation without steam in- 
jection. In this approach, a baseline condition is defined for a compo- 
nent that operate at its baseline conditions. For instance, in a gas tur- 
bine plant the baseline conditions are standard ambient The 
maintenance factor for baseline conditions is equal to one. At different 
operating conditions, maintenance factor is used to determine the in- 
creased level of required maintenance. Actually, the maintenance factor 
shows the impact of different operating conditions on system aging and 
consequently system maintenance intervals. The maintenance factor 
value is obtained from system historical data [21]. The sum of the 
hourly maintenance factors determines the actual age of a system that is 
known as EOHs [22]. Performance deterioration is a function of system 
actual age that can be derived from engineering experience and system 
historical data [23]. In fact, EOH concept helps to show the operating 
condition effect on system aging consequences such as performance 
deterioration and maintenance costs. 
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Fig. 7. Gas turbine flow diagram. 


Net output 314 MW 
Net heat rate 9422 kJ/kWh, LHV 


3. Problem formulation 
3.1. General description 


Scheduling is essentially a multidisciplinary task involving system 
performance, operations, maintenance, reliability, economics, and en- 
vironment dynamics [24]. Many of these factors must be considered in 
the modeling to have a comprehensive scheduling optimization model. 
This paper focuses on the performance, operations, and economic fac- 
tors that affect the power plant profit. The scheduling model optimizes 
operation and maintenance intervals of a power plant in order to 
achieve maximum power plant lifetime profit. In addition, components 
aging is a crucial factor in maximizing the profitability of power plants 
in long term [3], that is considered in the optimization procedure. In 
this section, the formulation of the proposed framework is described in 
detail. The section begins with the introduction of the model inputs and 
outputs. Then the formulation of the objective function is proposed, 
afterward the various constraint families are developed. 


3.2. Model inputs and outputs 


Three major external factors are the fuel cost, electricity price, and 
ambient conditions which are the dynamic factors over time and have 
impact on the final optimal solution. These parameters are functions of 
hourly, monthly and yearly time scales and are stochastic in nature. 
However, the trend of these parameters are important in long term 
scheduling that is deterministic. For instance, the price of electricity is 
higher in peak hours and this trend is deterministic in long term. The 
matrix of inputs for the year (k = 1), G) months and (i) hours is as Eq. 


(1): 


hx ft x2 h x12 
bx Lyx hx12 , , 

Iixjxı = . . L= ate JS eel ksl 
baxı b4 x2 Dax12 


(1) 


The model outputs are optimum operating conditions and main- 
tenance intervals which are derived from the optimization procedure. 


Oixı O1x2 O1x12 
O2x1 02x2 O2x12 

Oixjx1 = . : . b= l 22A FS 1212 ka 
Orax1 Orax2 Orax12 


(2) 
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Fig. 8. Variations of ambient temperature. 


Fig. 9. Variations of electricity price. 
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3.3. Objective function ie a 
Max Z= >) 2 (Rijk—Cijkm) 
The objective function maximizes the plant profit over the total a (3) 


operation life time horizon. The profit terms are system incomes (R), System incames ancl ade selling onnu power O ad aisefal eat 
and system costs (C). The profit for (k) years, (j) months, (i) hours, and (h) ss 


(m) maintenance types can be calculated as Eq. (3). 
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Gas turbine (GE 9F.05) maintenance intervals and costs. 


No. Type of inspection Interval [25] ERP [32]  Levelized aging cost (C%) 
- EOH % $/EOH 

1 Combustion 12,000 2.1 95.8 

2 Hot Gas Path 24,000 3.9 88.98 

3 Rotor 144,000 31.6 120.16 

k 2 24 

— P h 

R=), Pye X Ch + hir X Ch 
k=1 j=1 i=1 (4) 


(CP) and (C") are the electricity and heat unit prices, respectively, 
which change on hourly, monthly and yearly time scales. The system 
costs are the cost of fuel consumption and components aging cost. 


k 12 24 
1 j=l i=1 


k= 


C= (m x Ch + > EOHijkm x cs 


m=1 


(5) 


In Eq. (5), the first term is fuel cost and is calculated based on plant 
heat rate (hr) and unit fuel cost (Cf). The second term is aging cost and 
is defined as hourly preventive and corrective maintenance costs. (C“) is 
the cost according to the plant maintenance type (m) per its EOH. For 
instance, the (C*) of a gas turbine that has a rotor inspection cost of 
17,300,000 $ after 144,000 EOHs, would be 120.16 $/EOH. 

The EOHs, or equivalent operating hours, quantify the stresses that 
the system experiences asa result of generating energy and will in- 
crease according to variables including operating conditions, startups 


and peak operating hours and is calculated from Eq. (6), [22]. 
EOH = mf xt (6) 


(t) is hours of system operation and (mf) is maintenance factor. 
Assuming a day sample for each month, the EOH would be as Eq. (7). 
k 


12 24 
EOH = >) >) >) fine X 4; 
=1 i=1 


k=1 j 


(7) 


(dj) is the number of days in month (j). 


1s 
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Fig. 10. Variation of fuel cost. 
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Fig. 12. Validation of model results (solid lines) against experimental results (dotted 
lines) for aging effects on gas turbine output power and heat rate. 


3.4. Constraints 


The feasible region of the optimization problem is defined by dif- 
ferent constraints. Inequality constraints arise from limitations of 
system operations, and characteristics. An equality constraints are 
given by the equations which describe components performance func- 
tions and components aging. 


3.4.1. Operations and characteristics constraints 

Each system has an operating limit to be maintained. For instance, 
the constraint of an operating condition (0p;,) with a lower (op) and 
upper (OP) bands would be as Eq. (8). 

l 

Pig S Pjk S Pij (8) 

In addition, the preventive maintenance constraint should be ap- 
plied to the scheduling optimization problem. This constraint re- 
presents that the equivalent operating hours of the system must be 
lower than the system failure EOH which is presented in the standard 
maintenance policies. The maintenance constraint for the (k) years, (j) 
months, and (i) hours operation is as Eq. (9), [25]. 


Fig. 11. Validation of model results (solid lines) against 
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In the equation, (d;) is the day number of the month (j). 


3.4.2. Components I/O model 

The off-design behavior of each component is expressed by means of 
equations that relate its performance (outputs variables such as power 
generation and efficiency) to its operating conditions (inputs variables 
such as fuel and ambient conditions). An I/O model flow data is illu- 
strated in Fig. 3. 

The I/O model of a system is constituted by Eq. (10). 


P =f (A,F,S) (10) 


System performance (P) is a function of ambient conditions (A), fuel 
consumption (F) and system operation schedule (S). System perfor- 
mance can be the output power (p), efficiency (7), heat rate (hr), or 
maintenance factor (mf). 


P = [p,n,hr,mf | 


Because of the inputs variation, system I/O model outputs vary on 
hourly, monthly and yearly basis as well. 


(11) 


3.4.3. Components aging model 

In order to quantify component aging based on the system operation 
history, EOH approach is introduced. EOH is increased according to 
system operating conditions. In baseline condition the EOH of one-hour 
operation would be equal to one. For other operating conditions, EOH 
will be calculated based on the maintenance factor as Eq. (7). Plant 
performance deterioration (dP) is a function of EOH [23]. This 
function reflects the adjusted level of plant performance considering 
aging effects. 
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Fig. 13. Optimal operation schedules for different maintenance in- 
tervals. 
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dP% = f (EOH) (12) 


Fig. 4 represents the data flow of the components I/O and aging 
models integration. 

As is shown in Fig. 4, the aging based model represents deteriorated 
performance over time asa function of operating conditions. This 
function can be used in the scheduling optimization procedure in order 
to consider aging effects on optimal solution. 


4. Optimization algorithm 


The aging based optimal scheduling models are usually complex 
models which require a high structured solution methodology. 
Generally long term dynamic optimization methodologies are preferred 
in order to reach a deeper understanding of the process of aging in long 
term and determining optimum schedule [26]. In this study, a mixed 
integer nonlinear programming (MINLP) model which includes both 
continuous and binary variables is employed to estimate plant optimal 
schedule. The solution method is based on a hybrid heuristic optimi- 
zation algorithm of BARON and gradient search method. The branch 
and reduce optimization navigator (BARON) global optimization 
strategy is used to optimize plant continuous operation schedule. This 
methodology integrates conventional branch and bound with a wide 
variety of range reduction tests which most of them are based on 
duality [27]. The branching algorithm of BARON method is presented 
in Fig. 5. A random set of partition elements (open nodes) are gener- 
ated. One of the open nodes is selected and its related relaxation is 
solved. If the lower bound shows that the node is inferior than the 
current best known solution, the node is deleted and another node is 
chosen. If the node cannot be abounded immediately, it is divided into a 
set of new nodes that replace it on the list of open nodes; this process is 
branching. The procedure of node selection, bounding and branching is 
repeated until all the open nodes are abounded. 

At this step, the feasible region is branched into several smaller 
feasible spaces in which the model will be solved. The reduced feasible 
spaces help to improve the algorithm solution time and accuracy. The 
model is then solved for all feasible spaces and the global optimal so- 
lution will be derived. 

The output of the BARON optimization algorithm is long term op- 
eration schedule as is shown in Fig. 6. The next step of optimization 
algorithm is based on search method. In this stage, based on the system 
maintenance interval and maintenance cost, the best outage plan is 
derived to maximize lifetime profit of the system. Then results are 
checked if they meet the stopping criteria. Finally, the optimum deci- 
sion variable elements including plant operating conditions and outage 
plan will be derived. 


5. Case study 
5.1. Plant description 


The proposed framework is employed to determine the optimal 
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Fig. 15. Optimal hourly schedule for a typical day of each month of the first year of operation. 


schedule of a gas turbine power plant that serves a stationary load. The 
optimal schedule includes optimal operating conditions (hourly load 
and steam injection) and maintenance intervals during 15 years of 
operation. 

The studied gas turbine is a 9F.05 model of GE products. The flow 
diagram of the system is presented in Fig. 7. 

In addition, the main characteristics of the studied gas turbine are 
presented in Table 1 and more details of the system are provided in Ref. 
[28]. 
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5.2. Input profiles 


The hourly variations of the ambient temperature for 12 types of 
days are gathered from Tehran/Mehrabad weather station and are 
shown in Fig. 8, [29]. As it is expected, the ambient temperature is 
relatively low in the early morning, keeps increasing until noon, then 
decreases and reaches its minimum at midnight. The monthly varia- 
tions show that the average ambient temperature is relatively low in 
January, keeps increasing till August, and then decreases and reaches 
its minimum again in January. 

Fig. 9 presents the variations of the electricity price [30]. As can be 
seen, the price of electricity is higher at peak-loads than off-loads hours. 
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Fig. 16. Optimal hourly schedule for a typical day of each month of the second year of operation. 


The monthly variation of fuel cost is presented in Fig. 10, [31]. The 
figure shows that in cold months, natural gas price is lower than warm 
months which is due to higher fuel demand for heating in cold months. 


5.3. Scheduling formulation 


The aging based optimal scheduling for profit maximization, relies 
on the knowledge about economic performance of power plants. A 
projection of the effect of operating conditions and maintenance in- 
tervals on power plant output and performance deterioration are ne- 
cessary to evaluate power plant costs and incomes and consequently 
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power plant profit [15]. The objective of the optimization as presented 
in Eq. (13) is to maximize the long term profitability (Z) of a gas turbine 
power plant. The decision variables are operating conditions and 
maintenance intervals which are varying on hourly (i = 1,2,3....,24), 
monthly (j = 1,2,3.,...,12) and yearly (k = 1,2,3.,...,15) basis . 


15 12 24 
Max Z= > > » (Rijx—Cix ) 


k=1 j=1 i=1 


(13) 


Gas turbine plant income (Rj) is the income of selling electricity to 
the grid. 
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Fig. 17. Optimal hourly schedule for a typical day of each month of the third year of operation. 


(Pix) is the output power and (Cio) is the hourly unit price of elec- 
tricity. Moreover, plant costs (Cj) are the cost of fuel consumption and 
aging cost (maintenance hourly cost). 


3 


m=1 


(15) 


(hi;,) is the gas turbine heat rate and (Ci) is the monthly unit cost of 
natural gas. (mj;,,,) is the maintenance factor of inspection type (m). 
(C“) is the levelized cost of aging according to the inspection type (m). 
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Table 2. presents the inspection cost as a percentage of equipment re- 
placement price (ERP) for three types of inspection [25,32]. 

As can be seen, all terms of the profit are affected by components 
aging. Generated power and heat rate deteriorates over time due to 
components aging. Moreover, the aging cost reflects hourly cost of 
components aging. 

To show the effect of the operating conditions on profit terms dif- 
ferent coefficient factors are defined. The plant output power deviates 
from the base load output power (P°) of gas turbine because of the 
change in ambient temperature, steam injection, operating load and 
components aging as Eq. (16), [33]. 
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Py =P°xXaxBxaxd (16) 


In this equation, (@), (£), (6), and (A) represent the correction fac- 
tors of temperature, steam injection, system load percentage and aging, 
respectively. These correction factors can be calculated as Eqs. 
(17)-(19), [34-36]. 


a = —0.007 x T + 1.1076 (17) 
B = (1 + 0.042857 x %m,) (18) 
A = 0.03exp(—4.9 x 1074 x EOH) + 0.97exp(—5 x 1077 x EOH) (19) 


(T), (ms), and (EOH) are the ambient temperature, steam injection 
percentage, and equivalent operating hour. These functions are derived 
from system historical data [34-36]. 

In addition, plant heat rate deviates from its standard value (HR?) 
and can be calculated for different operating conditions based on the 
Eq. (20), [33]. 


hij = HR®°xa’x px’ xd (20) 


The (a@’), (6’), (6’), and (A’) represent the correction factors of tem- 
perature, steam injection, system load percentage, and aging, respec- 
tively. The correction factors are formulated as Eqs. (21)(24), [34-36]. 


a’ = 3 X 1075 x T? + 1.2 x 107°? x T + 0.9789 (21) 
B’ = (1—0.021429 x %m,) (22) 
6’ = —0.00026° + 0.003664—0.03226 + 0.13667—0.30236 + 1.365 (23) 
A’ = 1.02exp(3.4 x 10-7 x EOH)—0.02exp(—3.9 x 1074 x EOH) (24) 


Eqs. (21)-(24) are resulting from historical data of a gas turbine 
power plant [34-36]. In addition, GE company has developed main- 
tenance factor functions for heavy duty gas turbines. The functions 
determine maintenance factors based on system load and steam 
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Fig. 20. System hourly cash flow over calendar lifetime. 


Table 3 
Comparison between two optimal operating schedules. 





Aging based Base case Increase/decrease 
schedule schedule percentage (%) 
Profit (M$/yr.) 109.17 105.1 4 
Costs (M$/yr.) 55.14 55.88 -1 
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Fig. 21. System profit, income and costs for different interest rates. 


injection. The maintenance factors for combustion, hot gas path, and 
rotor inspections are calculated as Eqs. (25)-(30), [25]. 


mf! (combustion inspection) = k! x s (25) 


kt = Max (1,exp(0.34(%m;—%1))) (26) 


Fig. 19. Actual and calendar age of the studied gas turbine. 
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Fig. 22. Optimal hourly schedule for a typical day of each month of the first year of operation under 15% interest rate. 


1 up to base load 


-i 


exp(0.023 x (10.25 x %8—1025)) (27) 
mf? (hot gas path inspection) = k? x s (28) 
k? = 1 + (0.55 x %m,) (29) 
2, k 
mf? (rotor inspection interval) = 
1, non—peak (30) 


As shown, each maintenance factor depends on gas turbine 
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operating conditions including load (6) and steam injection (m,). The 
maintenance factor reflects the effect of operating conditions on system 
actual age (EOH) as Eq. (6). 

In addition, the case study model has several constraints. For in- 
stance, it is assumed that the gas turbine plant would operate until 
reaches its inspection intervals as presented in Table 2, [25]. 


m 


mf X di < EOH™ m= 1,2,3 


(31) 


Moreover, it is assumed that gas turbine startup takes 2 h and will 
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Fig. 23. Optimal hourly schedule for a typical day of each month of the second year of operation under 15% interest rate. 


consume more fuel until system reaches nominal operating status. 


Ôi—1 X 0i41 X E S Oj (32) 


where £ is a small positive infinitesimal quantity. This constraint shows 
that system cannot be offline just for one hour and the minimum offline 
hours in a row is two (startup time). 


6. Results and discussion 


The solution procedure is carried out based on a hybrid heuristic 
optimization algorithm. The general algebraic modelling system 
(GAMS) software is used to obtain the optimal values that maximizes 
gas turbine plant lifetime profit. The GAMS is a major software mod- 
elling facility that houses various software solver tools such as BARON 
(271. 
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In this section, results of the gas turbine model are presented and 
the effect of temperature, steam injection and EOH on system output 
power and heat rate generation is analyzed. Then, the studied gas 
turbine schedule (maintenance intervals and operating conditions in- 
cluding load and steam injection) is optimized for 15 years of system 
operation using proposed framework. Finally, performance deteriora- 
tion and economic evaluation of the studied gas turbine operating at the 
optimal schedule are presented. 


6.1. Validation 


In this section, the I/O and aging models of the gas turbine are 
validated. The output power and heat rate are considered as two cri- 
teria of gas turbine performance. Fig. 11(a) shows the correction factors 
of gas turbine output power (a) and heat rate (@’) asa function of 
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Fig. 24. Optimal hourly schedule for a typical day of each month of the third year of operation under 15% interest rate. 


ambient temperature. As ambient temperature increases, the com- 
pressor work increases and as a result output power decreases and heat 
rate growths. The model outputs are validated using experimental data 
gathered by [34]. As is shown in the figure, the model outputs provide a 
reasonable fit to the experimental data. 

Steam injection is a technique for gas turbine power augmentation. 
The effect of steam injection on output power and heat rate correction 
factors (8,6’) is illustrated in Fig. 11(b). as it is presented, seven percent 
steam injection can increase heat rate and output power by around 15 
and 30 percent, respectively. In the figure the comparison between 
model output and experimental results published by [35] are per- 
formed. Based on the results, the biggest difference observed is under 
3% that is negligible in our study. 

Performance deterioration as a function of EOH is used to validate 
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aging effects on performance deterioration over time. Fig. 12 presents 
the output power and heat rate correction factors (4,4) through system 
equivalent operating hours. As can be seen, the gas turbine output 
power drops and heat rate increases more quickly at low EOHs, as the 
gas turbine performance is very dependent on component aging and the 
degradation rate is higher at the beginning of system operating hours. 
In addition, as it can be seen, the model prediction is in good agreement 
with the experimental data published by [36]. At higher EOHs, model 
indicates lower output power and heat rate deterioration in comparison 
with experimental data. However, the modeling results are nearly close 
to the experimental data in this range as well. 
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6.2. Optimal maintenance scheduling results 


The optimum operating schedule depends on the continuous op- 
eration time span and maintenance interval. Fig. 13 indicates that op- 
timum operating conditions is changed due to different maintenance 
intervals. However, it is repeated for a specific maintenance interval 
during 15 years of operation. 

In order to find the optimum maintenance and operation schedule, 
the model is solved for 1-5 years of maintenance intervals and optimum 
operating schedule and plant indices including annual income, annual 
profit, annual fuel cost, and annual generated power are derived. The 
normalized expected indices as a function of maintenance intervals are 
illustrated in Fig. 14. It is found that as the maintenance inspection is 
postponed, the profit increases, until it reaches a maximum level. 
Afterward, the profit keeps decreasing if the maintenance inspections 
are postponed further. The results show that maximum profit is 
achieved when the maintenance is performed in the 3" year of system 
operation. 


6.3. Optimal operation scheduling results 


In this section, the optimized operation schedule for the derived 
optimum maintenance interval is presented. Figs. 15-17 show the op- 
timum hourly gas turbine load and steam injection for three years of 
operation. As can be seen, daily load and steam injection follow almost 
the same trend seasonally. As Figs. 9 and 10 show, the electricity price 
and natural gas cost profiles are the same for fall and winter. As a result, 
the optimum operation profiles in these seasons are almost the same as 
well. As it is shown, load is higher in the middle hours of the day. It 
increases as it goes from beginning hours into night and decreases in the 
last hours of the day. Steam injection is low at the first hours of the day 
and increases during peak hours and decreases in the last two hours. In 
spring and summer, the load and steam injection are low at the first 
hours of the day and increases gradually till midnight. This trend is due 
to the dynamics of electricity price and fuel cost in these seasons. 

As indicated by annual comparison, the power plant is scheduled to 
operate in a relative low output level with low power augmentation 
during the first months of operation. However, as shown, the gas tur- 
bine is turned on to the highest output level with steam injection at 
highest level during the last months of operation. 


6.4. Performance deterioration evaluations 


Performance deterioration through time is an important issue in 
evaluating gas turbine power plants operation. All power plants ex- 
perience losses in performance over time, even under normal operating 
conditions. However, performance losses can be optimized to have a 
lower economic impact. In the proposed framework, the performance 
deterioration is optimized in a most economical way. In this section 
performance deterioration for the derived optimal schedule (operation 
and maintenance) is presented. Fig. 18 depicts the output power and 
heat rate deterioration during system operating lifetime. As shown in 
Fig. 18, the performance deteriorates rapidly at the beginning hours of 
operation due to higher degradation rate in the first hours of operation 
[37]. Afterward the deterioration rate increases gradually that is be- 
cause of the smooth degradation rate and maintenance factor growth. 
Results indicate that after three years of operation, output power and 
heat rate deteriorate at about 10 and 7 percent, respectively. 

The relationship between gas turbine aging and calendar time is 
investigated in Fig. 19. The age of a gas turbine depends on factors 
including the chronological age, ambient conditions, number of startup, 
load, and steam injection. Determining the age of a gas turbine is fur- 
ther complicated by the fact that individual components age differently 
depending on these factors. However, the average age of a gas turbine 
can be evaluated based on the maintenance factor and is equal to the 
accumulative equivalent operating hours of the system as presented in 


1361 


Applied Energy 205 (2017) 1345-1363 


Eq. (7). As can be seen in Fig. 19, the gas turbine age growth up to 
140,000 EOHs through three years of operation. The actual and ca- 
lendar age difference increases as time passes. For instance, the actual 
age is around 3, 4, and 5 times higher than calendar age after one, two 
and three years of operation, respectively. In fact, the ratio of the actual 
age to the calendar age at time (t) shows the average hourly main- 
tenance factor of the gas turbine operating until time (t). 


6.5. Economic evaluations 


Fig. 20 depicts the cash flow of the gas turbine power plant. In order 
to analyze system profit, the income of system and costs are con- 
currently illustrated in Fig. 20. Profiles fluctuation is because of the 
seasonal variations of the market signals and optimal operation sche- 
dule. According to the results, the average income is low in the be- 
ginning and end periods of each year, due to the lower electricity price. 
At lower electricity price, system prefers to generate power at lower 
rate to reduce components aging and fuel costs. However, from the 
beginning of July to the end of September, the electricity price is higher 
and the income of generating power outweigh aging and fuel cost. 
During period April to Jun, the income is in its lowest levels if com- 
pared with other months because of reducing energy market price. 
Furthermore, it is found that the scheduled system cost follows the 
same trend as that of fuel cost in a year. It increases as it goes from 
winter into summer and decreases from fall to spring. Furthermore, the 
fuel cost profile is repeated annually and its average value increases 
over time as a result of components aging and heat rate deterioration. 


6.6. Aging consideration effectiveness 


The optimized aging based schedule is compared with the base case 
schedule (full load without steam injection) in terms of economical 
comparison that is summarized in Table 3. In fact, base case schedule is 
the output of a scheduling optimization model without considering 
aging effects. Results indicate that the aging based optimal schedule 
leads to a more efficient and economic operation of gas turbine than the 
base case schedule. Based on the results, during 15 years of operation, 
the operation costs are nearly the same. However, operating at aging 
based optimal schedule increases gas turbine annual profit up to 4% 
that shows the effectiveness of proposed framework. 


6.7. Sensitivity analysis on interest rate 


Since the problem of optimal scheduling is a long term optimization 
problem, it is more realistic to take the uncertainties of investment 
opportunities into consideration. In this study, three different interest 
rates are considered in the optimization model. 

The optimal operating and maintenance schedule is derived for 
three different interest rates of 5, 10 and 15%. Fig. 21 shows the nor- 
malized income, cost, and profit of the plant under different levels of 
interest rates. As can be seen, the net present values of profit, income, 
and costs decreases by increasing the interest rate. 

Figs. 22-24 show the optimum hourly load and steam injection for 
three years of operation under 15% interest rate. It should be noted 
that, the model is solved for 1-5 years of maintenance intervals and 
optimum operating schedules are derived. The results show that max- 
imum profit is achieved when the maintenance is performed in the 3" 
year of system operation. 

Results show that considering interest rate changes optimal sche- 
dule entirely. Increasing interest rate leads to an operation schedule 
with higher steam injection and lower load at the beginning hours of 
operation. The increase in steam injection level and load are the ways of 
power augmentation. Steam injection can increase power with lower 
fuel consumption in comparison with load increment. However, has a 
higher degradation rate and deteriorates system performance more 
rapidly. 
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In the case of ignoring interest rate, system operates at the operation 
schedule with lower degradation rate over time. For instance, the 
power is increased at the first months of operation. Conversely, the 
power augmentation is handled by steam injection at the end of the life. 

On the other hand, under high interest rates, system prefers to op- 
erate at higher income and lower fuel cost at the beginning hours of 
operation despite the increase in performance deterioration over time. 
It is because of the lower value of cash flow over time under high in- 
terest rate environment. 


7. Conclusions 


In this paper an aging based optimal scheduling framework is pro- 
posed to maximize power plants lifetime profit. The framework con- 
siders components aging effects in the optimization procedure. The 
equivalent operating hour approach is introduced to determine the ef- 
fect of component aging on system performance deterioration, aging 
cost and consequently optimal schedule. The aging cost is the hourly 
maintenance cost that is considered as a term in the objective function. 
This shows the dependency of maintenance cost on operation schedule 
using maintenance factor coefficient. 

The framework is developed based on the mixed integer nonlinear 
dynamic optimization approach. A two stage optimization algorithm is 
used as the solution methodology. The first stage is based on the 
BARON method and determines the continuous optimum operating 
schedule and the second stage is a search method that provides the 
optimal maintenance schedule. 

Finally, the proposed methodology was applied to a gas turbine to 
evaluate the framework effectiveness. Results show that the EOH ap- 
proach is a promising method to predict gas turbine performance as it 
degrades over time. Moreover, the obtained results from scheduling 
optimization model show that the aging based scheduling framework 
can be an appropriate approach to achieve maximum lifetime profit. 
Through the presented analysis, it is shown that operation scheduling 
and maintenance intervals should be optimized simultaneously. 
Furthermore, when compared to the base case schedule (optimal 
schedule without considering aging effects), the aging based framework 
is able to achieve higher profit up to 4% annually. 

Furthermore, a sensitivity analysis on interest rate is performed. It is 
shown that optimal schedule depends on interest rate level strongly and 
in schedule optimization problems economic knowledge and assump- 
tions are essential. 

In this study, the framework objective function is power plant life- 
time profit. However, any other objective function, such as primary 
energy consumption, plant efficiency, or pollutant emissions, could be 
implemented in the framework here defined. 

The major concluding remarks can be summarized as follows: 


e Results demonstrate gas turbine power and heat rate deteriorate as 
time pass, and the economic impact of these issues are significant. 
Therefore, aging based optimal scheduling should be performed to 
achieve maximum lifetime profit. 

The optimization results show that the optimal operation schedule 
depends on the maintenance interval. Therefore, operation and 
maintenance schedule should be optimized simultaneously for the 
entire lifetime of the plant. 

To show the effectiveness of the developed framework, model out- 
puts are compared with an optimal case that doesn’t consider aging 
effects in the optimization procedure. Results illustrate operating at 
aging based optimal schedule is more beneficial for long term profit 
maximization and the power plant will have 4% higher annual 
profit. 
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